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ABSTRACT 



Context. The discovery of about 700 extrasolar planets, so far, has lead to the first statistics concerning extrasolar planets. The pres- 
ence of giant planets seems to depend on stellar metallicity and mass. For example, they are more frequent around metal-rich stars, 
with an exponential increase in planet occurrence rates with metallicity. 

Aims. We analyzed two samples of metal-poor stars (-2.0 < [Fe/H] < 0.0) to see if giant planets are indeed rare around these ob- 
jects. Radial velocity datasets were obtained with two different spectrographs (HARPS and HIRES). Detection limits for these data, 
expressed in minimum planetary mass and period, are calculated. These produce trustworthy numbers for the planet frequency. 
Methods. A general Lomb Scargle (GLS) periodogram analysis was used together with a bootstrapping method to produce the detec- 
tion limits. Planet frequencies were calculated based on a binomial distribution function within metallicity bins. 
Results. Almost all hot Jupiters and most giant planets should have been found in these data. Hot Jupiters around metal-poor stars 
have a frequency lower than 1.0% at one sigma. Giant planets with periods up to 1800 days, however, have a higher frequency of 
fp = 2.63^Q g%. Taking into account the different metallicities of the stars, we show that giant planets appear to be very frequent 
(fp = 4.48^|3g%) around stars with [Fe/H] > -0.7, while they are rare around stars with [Fe/H] < -0.7 (< 2.36% at one sigma). 
Conclusions. Giant planet frequency is indeed a strong function of metallicity, even in the low-metallicity tail. However, the frequen- 
cies are most likely higher than previously thought. 

Key words. Planetary systems - Techniques: radial velocities 



1. Introduction 

Since t he discove ry of the first extrasolar planet in 1995 (51 
Peg b, iMavor & O ueloz 1995), the search for extrasolar plan- 
etary systems accelerated. Today, around 750 planets are an- 
nounced. IVlost of them were detected using the radial veloc- 
ity technique. Although 750 is a relatively high number, the 
theory of planet formation and evolution is still under debate 
dPollack et al.lll996l: iMaver et al.ll2002l: iMordasini et al.ll2009h . 
The situation is particularly difficult for giant planet formation. 
Currently, there are two proposed models: core-accretion (e.g . 
iPollacket al.ll 19961: iRice & Armitaggi2003l: lAlibert et al.ll2004i) . 



where gas from the protoplanetary disk is accreted around a pre- 
vioti sly formed r ocky/icy core, and the disk instability model 
(e.g. lBosslll997l : iMaver et al]|2002h . where a planet is formed 



* The data presented herein are based on observations collected at 
the La Silla Parana Observatory, ESO (Chile) with the HARPS spec- 
trograph at the 3.6-m telescope (ESO runs ID 72.C-0488, 082.C-0212, 
and 085.C-0063) and at the W.M. Keck Observatory that is operated as 
a scientific partnership among the California Institute of Technology, 
the University of California and the National Aeronautics and Space 
Administration. This Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation. 
** Table 1 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qc at ?J/A-l-A/| 



because of a direct gravitational instability in the protoplane- 
tary disk, in the same way as stars form from inter stellar clouds. 
A help ful overview of l^oth models is given by iMatsuo et al.l 
(I2OO7I) . 

One of the main advantages of the instability model is the 
timescale that is needed to form the planets. Early results sug- 
gested that the slow accretion phase (about lOMyr) in the core- 
accre tion model may tak e longer than the lifetime of a T Tauri 
disk (IPollacket alJll996h . In that sense, giant planets could not 
form within the core-accretion model. However, more recent re- 
sults suggest that this may not be a real problem. In fact, it has 
been shown that the process can, for exam ple, be accelerated 
by including disk-indu ced orbital migration (lAlibert et al.ll2004t 
IMordasini et al]|2009l) . 

Theories of migration became more important with the 
discovery of 51 Peg b and other hot Jupiters. These close- 
in giant planets are highly unlikely to have formed in-situ. 
Interestingly, however, disk-induced migration does not neces- 
sarily provide the correct explanation for al l these hot Jupiters . 
Other ideas have been put forward (e .^ 



Triaud et 



t Ju piters . 
aP boiOl 



ira (e .g. 

iMorton & Johnsoiil 1201 ll ISocrates et al.l 120121) . which include 
more "violent" migration mechanisms. Discoveiies o f giant 
plane t s on wide orbits of tens to hundreds AU (Maro is et al.l 
I200I iLagrange et al.ll2010i) also raise questions. Overall there 
is still space to support that the disk-instability model may 
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be at work, at least to explain pa rt of the detected planets 
(IVorobvov & Basull2010H Boss"201 IV 

Additional clues about this problem come from the analy- 
sis of planet-host stars. The presence of a planet seems to de- 
pen d on several stellar pr operties, such as mass and metallic- 
ity (lUdrv & Santosll2b07h . Concerning metallicity, it has been 
well-established that more metal-rich stars have a higher prob- 
ability of harboring a giant planet than their lower metal- 
licity counterpart s (Gonzalez 1997; Santos et al. 2001 . l2004t 
Fisch er & Va lentil j2005L lUdrv & Santosi i2007i: i Sozzetti et al.l 



2009 ; Sousa etaLll201 la) . The occurrence rate even increases 



dramatically with increasing metallicity. Current numbers, based 
on the CORALIE and HARPS samples, suggest that around 
25% of the stars with twice the metal content of our Sun are 
orbited by a giant planet. This number decreases t o ~ 5% fo r 
solar-metallicity objects (Sousa et al. 201 lb; Mavor et al.ll201 ih . 
A similar trend was also obtained by previous results (e.g . 



ISantos et al.l2004l:lFischer & Valentil200llJohnson et al.l2010l) . 
Cu riously, no such trend is observed for the lower mass plan- 
ets (lUdrv et al.ll2006tlSousa et aLll200i iMavor et al.ll20TTl) . The 
Neptune-mass planets f ound so far seem to have a rather flat 
metal licity distribution dSousa et al.ll2008l 1201 Ibt iMavor et al] 
I20T1I) . 

This observed metallicity correlation favors the core- 
accretio n model f or the formati on of giant planet s (llda & LinI 
[2004; Udrv & Santo s 2007; Mordasin i etaP |20 1 2l) because the 
higher the grain content of the disk, the easier it is to build the 
cores that will later accrete gas. According to the disk-instability 
model, the presence of planets w ould not be strongly dependent 
on stellar metallicity (Boss 2002). 

Understanding the frequency of different types of planets 
around stars of different mass and metallicity is thus providing 
clues about the processes of planet formation and evolution. This 
has inspired the construction of specific sam ples t o search for 
planets around different types of sta rs (e.g. ISantos et al.il2007t 
ISato et al.ll2008l; ISozzetti et al.ll2009t) . Statistics of these samples 
will help in understanding the formation processes and constrain 
the models. 

We present an analysis of two metal-poor samples that were 
designed for planet-finding purposes. In Section |2] an overview 
is given of the samples and their data. Section lSTI reports on the 
detection limits of these samples. Planet frequencies are calcu- 
lated in Section |4] Finally, conclusions are made in Section |5] 
together with a discussion. 




I 1 HARPS 

cr: KECK 




stellar mass (M sun) 



Metallicity [Fe/H] 



Fig. 1. Relative histogram of the number of measurements (top 
left panel), rms noise (top right panel), stellar mass (bottom left 
panel) and metallicity (bottom right panel) in the two datasets. 
The blue line represents the HARPS sample, the green line the 
KECK-HIRES sample. Both graphs in the top panel are cut at 
20 for better visibility. The top left panel has 19 stars from the 
HARPS sample higher then 20. The top right panel has 1 star 
from the KECK sample higher then 20. 



high rotation). Most of the stars in the final sample of 88 stars 
have five or more measurements with an rms of ~ 1 - 2.5 m s"' 
as shown in Fig. [1] and Table [1] In the bottom panels of Fig.[T] 
the mass and meta llicity distribution is shown. Values were taken 
fromlS ousa et al.l f201 l a). The (spectroscopic) metallicities dif- 
fer slightly from the photometric estimate that was initially used, 
including a few outliers. All values are listed in Table [T] 

Three planetary mass companions were found in this sample . 
They are orbiting HD18172CI HD1 90984 (iSantos et al.ll2010l) . 
and HD171028 (Santo s et al]l2007h . AU three are giant plan- 
ets in long-period orbits. A f ourth planet candidate, orbiting 
HD107094, was announced in ISantos etaH (1201 ll) . With a 4.5 
Mjup minimum mass and a 1870 day period, it would again be 
a giant planet in a long-period orbit. However, this planet could 
not be fully confirmed yet. 



2. Data 

Radial velocity measurements from two different samples of 
metal-poor, solar-type stars were used in this paper. 

2.1. The HARPS sample 

ISantos et al.l (1201 ll) reported on the first sample. They observed 
104 metal-poor or m ild metal-poor sol ar-type stars with the 
HARPS spectrograph (Mav or et alll2003l) . The objects were ob- 
served from O ctober 2003 till July 2010. Based on the catalog of 
iNordstrom et al.l (l2004i) . aU late-F, G and K stars south of H-10° 
declination with a V magnitude brighter than 12 were chosen. 
After discarding spectroscopic binaries, giant stars and active 
stars, Santos and coworkers only recovered the 104 stars with 
an estimated photometric [Fe/H] between -0.5 and -1.5. After 
the observations, 16 stars were also discarded because they were 
unsuitable targets for planet-finding purposes (binarity, activity. 



2.2. The KECK-HIRES sample 

iSozzetti et al.l (l2009l) reported on the second sample. They ob- 
served 160 metal-poor solar-type stars with the HIRES spec- 
trograph on the Keck 1 telescope at Mauna Kea in Hawaii 
(IVogt et al.lll994l) . The objects were all observed at least twice 
over a timespan of three years (2003 - 2006). This sample of 
stars was drawn from the Carney -L atham and Ryan samples 
(e.g. ICarnev et al.llT994t lRvanlll989h . Additional criteria were 
applied, and in the end, Sozzetti and collaborators chose the 160 
stars with a V magnitude brighter than 12, an effective temper- 
ature Teff lower than 6000K and a metallicity [Fe/H] between 
-0.4 and -1.8. All stars are situated north of -25° declination. 
Most stars in this sample have 4 to 10 measurements with a rms 
of ~ 9 ms ' as seen in Fig. [1] and Table [1] Photon noise is the 
main contributor to these higher rms values. The m ass and metal- 
licity values, seen in this figure, are taken from ISozzetti et al.l 
(I2OO9I) . AU values are listed in Table [1] Typical uncertainties on 
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Table 2. Names and number of measurements of stars that are 
present in the two samples. 



name HARPS 




name KECK 




BDp062932 


4 


G66-22 


4 


BDp083095 


3 


G16-13 


5 


HD 104800 


6 


Gll-36 


4 


HD111515 


5 


G14-5 


7 


HD 126681 


13 


HD 126681 


6 


HD131653 


4 


G151-10 


7 


HD 134440 


10 


HD 134440 


6 


HD 1482 11 


31 


HD148211 


3 


HD148816 


6 


HD148816 


10 


HD193901 


3 


HD 193901 


5 


HD 196892 


3 


HD 196892 


5 


HD215257 


37 


G27-44 


8 


HD22879 


36 


G80-15 


5 


HD88725 


22 


G44-6 


5 



1.0 



0.8 




10° lo' 10^ 10^ 

Period (days) 



Teff , [Fe/H], and M» are 100 K, 0.1 dex, and 0.1 Mq respec- 
tively. 

No planetary signals were found in this sample. 
2.3. The combined sample 

Fourteen (14) stars have measurements in both samples. This 
amounts to a complete sample of 234 metal-poor solar-type 
stars. The two samples use different naming for the stars (see 
Table|2]). In this paper, the naming from the HARPS sample will 
be used for these 14 stars. 

Combining the measurements of these stars provides more 
data to look for possible planetary signals. The data of two dif- 
ferent telescopes were combined by subtracting the mean of the 
data from each set in the overlapping time-interval. A general 
Lomb-Scargle (GLS) peiiodogram analysis (see Section I3.lt 
was then performed on the 14 stars. No significant peaks were 
found in these periodograms. 

3. Detection limits 

3.1. Methodology 

In the literature, different authors used two main approaches to 
find detection limits in RV data. One is based o n y^- and F-tests 
(e.g. lLagrange et al.ll2009l:ISozzetti et al.ll2009l). another is based 
on a periodogram analysis (e .g.lCumming et al.ll999tlCumming| 
[2004; Endl et al. 200 li: iNaravan et al. 20051). In this paper, the 
second approach was chosen because we consider that we have 
enough measurements for a reliable periodogram analysis (see 
below for a comparison of the methods). 

A frequency analysis of unevenly sampled data (like RV 
data) can be performed by using the GLS periodogram (Scargle 
ll982LIZechmeister & Ki irster 2009). The GLS is equivalent to a 
least-squares fitting of a full sine-wave, including weights and an 
offset, representing a circular orbit. In the resulting periodogram, 
the power p{a>) is calculated as a function of frequency. This 
power measures how much the fit to the measurements improves 
by using a sinusoid instead of a constant. This analysis can also 
be performed by replacing the "sine-wave" with a Keplerian 
function. An example of a GLS periodogram for a Keplerian fit 
is shown in Fig. |2] (top panel). This plot refers to a GLS of the 
RV data of HD 134440, a star that has 16 measurements with an 
average rms of 4.9 ms"' and a timecoverage of 1885 days. We 
used both the circular and Keplerian approach. 




lO' 10^ 

Period (days) 

Fig. 2. Example of the GLS periodogram (top panel) and the de- 
tection limits (bottom panel) for the star HD 134440 with mea- 
surements from HARPS and KECK-HIRES. For both plots, a 
Keplerian fit was made. In the top panel, the power is plotted 
against the period. The horizontal solid line marks the power 
level for an FAP of 1%. The bottom panel plots the minimum 
planetary mass against the period. The solid line represents the 
detection limits for these data. The dashed lines indicate a circu- 
lar planetary signal with an RV semi-amplitude of 1, 3 and 9 m 
s"' (lower to higher line). 



The significance of a peak in the GLS periodogram can be 
determined analytically. However, not all values, like the num- 
ber of independent frequencies, are properly determined in this 
way (Zechmeister & Kiirster 2009). Alternatively, the signifi- 
cance can be ac cessed if we use a bootstrapping method (e.g. 
lEndl et al]|200II:lDumusque et al...201 1.) . Multiple time series of 
radial velocities are made by shuffling (with repetition) the real 
radial velocities while preserving the original times. On each vir- 
tual time series, a GLS is performed to determine the highest 
peak (frequency independent) in the periodogram. This can be 
used to determine the percentage of bootstrapped periodograms 
with maximum peaks above the one observed in the GLS of the 
actual data. This procedure allows us then to derive the false 
alarm probability (FAP) level (see top panel Fig. For this 
work, we have chosen to adopt 1000 bootstrapping series to es- 
timate the significance of the peaks. 
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Table 1. Relevant values for the targets in the two samples. The complete table is provided in electronic form only. 



Star 


n 


Mean rms 


Timespan 


Tf// 


[Fe/H] 


M. 




Reference 






[ms"'] 


[days] 


[K] 




[Mo] 






HD123517 


9 


3.02 


1596 


6082±29 


0.09±0.02 


1.214 


:0.08 


(1) 


HD 124785 


17 


2.01 


1518 


5867±21 


-0.56±0.01 


0.874 


:0.02 


(1) 


HD 126681 


13 


2.23 


1964 


5570+34 


-1.15±0.03 


0.714 


:0.02 


(1) 


G15-7 


6 


9.27 


890 


5280 


-0.88 


0.74 




(2) 


G151-10 


7 


10.91 


891 


5287 


-0.70 


0.76 




(2) 


G157-93 


3 


4.43 


115 


5409 


-0.78 


0.78 




(2) 



References. (D lSousa et all ( 1201 lah : (2^ ISozzetti et al.1 ( 120091) . 



Detection limits in RV data are derived by inserting a fake 
planetary signal in the data (circular or Keplerian). The proce- 
dure goes as follows. Virtual time series are made by adding 
these signals to the original data, which are treated as random 
noise. For a circular orbit, a fake signal 



y(t) = K sin 



In 

P ^ 



+ c 



(1) 



is added to the original data. Virtual series were made for 
periods P from 0.5 to 3000 days, semi-amplitudes K from to 
10 km s"' and ten phases ip, evenly separated by tt/ 10. On each 
series, a GLS periodogram is performed. For each period, a sig- 
nal is considered detected if the periodogram gives a peak at that 
period with a FAP of 1% for all 10 phases. The minimum semi- 
amplitude K for which a signal is detected expresses the lower 
limit for detectable planets in these data. 

The same approach can be taken for eccentric Keplerian sig- 
nals. The fake signal, added to the original data, can in this case 
be described as follows: 



y{t) - a cos v(f) -i- b sin v(f) -i- c 



(2) 



with a - K cos m, b — -K sin m and c - Ke cos m + y. Here, 
K is the RV amplitude, e the eccentricity, m the longitude of 
periastron, y the constant system RV and v(f) the true anomaly. 
This true anomaly is a function of t, e, P and Tq, the time of 
periastron passage. Again, virtual series were made for periods 
P from 0.5 to 3000 days and semi-amplitudes K from to 10 km 
s"' were tried. For each period and semi-amplitude, 1000 virtual 
signals were created, with 10 different eccentricities e (between 
and 1), 10 different longitudes of periastron zn (between and 
2n) and 10 different times of periastron Tq (between and P), all 
evenly separated. In this case, a planet at a specific period P is 
considered detected if the periodogram gives a peak with a FAP 
of 1% for all 1000 signals. 

The minimum semi-amplitudes can then be transformed into 
planetary masses (expressed in Earth mass) with the following 
formula: 



94 I ^TlPM 

M„ sin i = 7.4 ■ 10'^*K Vl - 1 



2\l/3 



(3) 



where the semi-amplitude K is expressed in m/s, the period 
P in days, the stellar mass M, in kg and the gravitational con- 
stant G in m-'kg^'s"^. An example of these detection limits for 
HD 134440 is shown in Fig. |2](bottom panel). A Keplerian signal 
was inserted to obtain these limits. 

As a comparison, for some stars, the detection lim- 
its were also calculated following the method described in 



Lagrange et a(. (2009) method 
GLS method 




10 

Period (days) 



Fig. 3. Planetary mass is plotted against period for G 19-27 (Keck 
sample). The blue (upper solid) line shows the detection limits 
for circular planetary signals in the data with a FAP of 1%. The 
black (lower dashed ) line shows the detecti on limits based on the 
method described in lLagrange et al.l (l2009h . 



iLagran^e et alJ (l2009l) . Virtual RV sets are created with the ex- 
pected RVs for a circular orbit, added with a random noise be- 
tween + RV error, where the RV error is the mean error of 
the real data. For every given period considered (between 0.5 
and 1000 days), 200 virtual data sets were taken by varying the 
phase of the signal. This was then performed for different semi- 
amplitudes till the signal was detected. A signal was considered 
detected if the standard deviation of the real RV measurements 
was less than the average value of the virtual standard deviations. 
This resulted in detection limits with overall the same shape, but 
a factor of 2.5 lower, as can be seen in Fig. |3] This shift can be 
explained by the fact that the periodogram analysis used in this 
work is more conservative in its definition of detectability, be- 
cause it needs a peak above the 1% FAP level for all phases. 
Because the overall shape is the same, the conservative peri- 
odogram analysis is favored for the purposes of this work. 

Note that in all these analyses, a circular (resp. Keplerian) fit 
to RV data needs at least four (resp. seven) measurements. To 
be more conservative, the choice was made for at least six (resp. 
ten) measurements for the analysis. 
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HARPS sample 




10 10 
Period (days) 

KECK sample 



E lo' 




lO' 10' 
Period (days) 

Combined sample 




10 10 
Period (days) 



Fig. 4. Planetary mass is plotted against period. The blue line 
shows the detection limits for circular planetary signals in the 
data with a FAP of 1%. In the top panel, the data from the 
HARPS sample are shown. In the middle panel the data from 
the KECK-HIRES sample are shown and in the bottom panel, 
the combined sample is shown. 



3.2. Stars with at least six measurements 

Detection limits based on a circular fit were made only for stars 
with at least six measurements. This resulted in 64, 50, and 1 14 
stars (72.7, 31.25, and 48.7%) for the HARPS, KECK-HIRES, 
and the combined sample, respectively. 

The long-period signals (longer than the timespan of the 
measurements) produce a high-amplitude power in the GLS pe- 
riodograms. To analyze the data for the presence of "shorter" 
period peaks, it is consequently necessary to remove these. A 
linear function was fitted to the RV data of each star, using a 
least-squares fit. If the correlation coefficient r was greater than 
0.7, the linear fit was considered relevant and subtracted from 
the original data. This was the case for HD107094, HD11397, 
HD215257, HD123517, HD88725, HD144589, and HDl 13679 
in the HARPS sample, G135-46, G63-5, G197-45, HD134439, 
G237-84, HD192718, HD215257, HD7424, and G63-44 in 
the KECK-HIRES sample and all these objects together with 
HD193901 in the combined sample. The planetary signals from 
the three confirmed planets in the HAR PS sample (with the pa- 
rameters taken from lSantos et al.l201 lb were also subtracted, be- 
cause we did not aim to confirm their existence but rather to 
search for additional signals. 

Figure |4] displays the detection limits for circular planetary 
signals in the three samples. Minimum planetary mass is plotted 
against period. The limits shown in the figure correspond to the 
values for which a signal can be detected in 80% and 95% (blue 
and green curve, respectively) of the stars in the sample. For 
clarity, we also include two dashed lines at Jupiter mass and at 
50 Me, a value close to the generally accepted lower limit for 
giant planets. From these plots, it is clear that at a 80% level 
no hot Jupiters (here defined as having an orbital period <10 
days) could have been missed in the HARPS sample, though 
none was detected. In the KECK and the combined sample, most 
hot Jupiters should have been found. 

3.3. Stars witli at least ten measurements 

For stars with at least ten measurements, the detection limits 
were also calculated for Keplerian signals. This was mainly use- 
ful for the HARPS sample, where more stars have dozens of 
measurements. In the HARPS, KECK-HIRES, and combined 
sample, respectively, there are 37, 7 and 47 stars (42, 4.4 and 
20.1%) with at least ten measurements. 

Figure |5] displays the Keplerian detection limits for the 
HARPS sample. The same trends and planetary signals as for 
the circular fit were subtracted before the analysis. The limits 
are for a 80% and 95% sample completeness (blue and green 
curve, respectively). The limits for the KECK sample are not 
shown because there are too few stars in the sample with at least 
ten measurements. Again, it is clear that no hot Jupiters were 
missed in the HARPS sample. The same is true for the com- 
bined sample. Note also that our data in the HARPS sample are 
sensitive to the detection of planets with masses above that of 
Jupiter for the whole covered period range. 



4. Planet frequency 

With these limits and the number of planets found in the sam- 
ples, a statistical analysis can be made of the giant planet fre- 
quency as a function of stellar metallicity. There are again sev- 
eral approaches to perform this analysis, among which the two 
main ones are a binning approach and a parametric approach 
(see below). The former makes use of a binomial distribution. 
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50 Me,,,,, 







10° 10^ 10^ 10^ 

Period (days) 



Combined sample 



10^ 























10° lo' 


10^ 


10^ 



Period (days) 

Fig. 5. Planetary mass is plotted against period. The blue line 
shows the detection limits for Keplerian (non-circular) planetary 
signals in the data with a FAP of 1%. In the top panel, the data 
from the HARPS sample are shown and in the bottom panel, the 
combined sample is shown. 



The probability of finding n detections in a sample of size can 
be calculated as a function of the true planet frequency fp-. 

Th is method is described in the appendix of iBurgasser et al.l 
(1200 3^. For this asymmetric distribution, the errorbars can be 
computed by measuring the range in fp that covers 68% of the 
integrated probability function. This is equivalent to the 1-sigma 
errorbars for a Gaussian distribution. 

As seen above, no hot Jupiters were found in our samples, 
while the detection limits indicate that they most likely should 
have been detected. Zero (0) detections in a sample of 1 14 stars 
(stars in the combined sample with at least six measurements) 
leads to a frequency /p < 1.00% (calculated from/;, = 0.37+j}^). 
The frequency rises to fp < 2.36% if only the 47 stars with at 
least ten measurements are taken into account. 

For giant planets in general (i.e. planets with a mass higher 
than SQMEartii), there are three detections in a sample of 114 
stars, which gives a frequency f, - 2.63^^ g% (see Fig. |6ll. In 
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Fig. 6. Probability as a function of true planet frequency for a 
given amount of detections n and sample size A^. The solid verti- 
cal line denotes the observed planet frequency, while the dashed 
lines show the limits of the centered 68% area, thus expressing 
the 1-sigma errorbars. 



this sample of stars, 90% have timespans longer than 900 days. 
This makes these frequencies sensitive to planets with periods up 
to 1800 days. If the sample is limited to only the 47 stars with at 
least ten measurements, the frequency becomes f, = 6.38^^ q%. 
For this smaller sample, this frequency is sensitive to planets 
with periods up to 2600 days. 

Our sample has a metallicity distribution that peaks around 
-0.7 dex. If we divide it into two parts (above and below this 
limit), we find that on the high-metallicity side there are 67 (resp. 
35) stars with at least six (resp. ten) measurements. This is also 
the side of the distribution of the detected giant planets. This thus 
leads to percentages off, = 4.48!f°^% (resp. = 8.57+^21%). 
Around the stars in our sample with metallicities lower than -0.7 
(47 stars with at least six measurements), no planets were de- 
tected. This again gives a frequency of fp < 2.36%. To check if 
the values are dependent on the choice of the bins, we repeated 
the calculation changing the position of the bins by 0.1 dex (the 
typical 3-sigma error of the individual metallicity estimates), 
thus changing the high (low) metallicity side of the distribution 
to cover the range of [Fe/H]> -0.8 dex (< -0.8 dex, respec- 
tively). In these bins, the planet frequencies for stars with at least 
six measurements become f, = 3.70;^j j3% resp. f, < 3.32%. 
Both results are comparable within the errorbars. 

Alternative ly, a parametric approach, similar to the one used 
bv lJohnson et al. (2010), was also considered. Here, the data are 
fitted with a functional form, dependent on stellar metallicity. As 
a functional form, we chose 

fp^C- Kf^'"'/"^ , (5) 

which is typically used for solar neighbourhood samples. 
The best parameters (fi, C) are then determined by using a nu- 
merical fitting procedure, based on Bayesian infe r ence. D etails 
of this procedure can be found in ^J ohnson et al.l (l2010l) . As a 
prior, the choice was made for uniformly distributed parameters 
over [0.0,3.0] and [0.01,0.30] for/? and C, respectively. The 
best parameters were found to be = 1.3 and C - 0.17. The 
mean metallicity for stars with six measurements and a metal- 
licity higher than -0.7 dex, is -0.55 dex. For the derived /3 and 
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C, the corresponding expected frequency would be 3.29%, com- 
patible with the results of the binning procedure. Owing to the 
limited size of the sample and because we are exploring a metal- 
licity regime not previously explored with high-precision radial 
velocities, it is unclear which appropriate functional form should 
be used. 

5. Conclusions and discussion 

Radial velocities of two samples of metal-poor solar-type stars, 
taken with two different instruments (HARPS and KECK- 
HIRES), were used to detect extrasolar planets. Only three gi- 
ant long-period planets were found out of the 234 stars, together 
with one giant candidate. Fourteen stars were present in both 
samples, but the expanded datasets and extended baseline of the 
observations did not reveal any additional signals. 

After subtracting linear trends or planetary signals, detection 
limits in the samples were calculated. The method was based 
on a GLS periodogram analysis and bootstrapping. Limits were 
calculated for circular and Keplerian signals. These lower limits, 
as shown in Figs. |4]and|5] are expressed in minimum planetary 
mass and period. For the stars in the KECK samp le, detection 
limits were already derived bv ^ Sozzetti et al.l (l2009h . They used 
a method based on;^'^- and F-tests. The detection limits derived 
in this work perfectly agree with their results. 

A statistical analysis was performed to estimate the planet 
frequency around metal-poor main-sequence stars. Taking into 
account only the stars with at least six measurements, we showed 
that the frequency of hot Jupiters around metal-poor stars is 
lower than 1.00%. T his is consistent with previous studies 
(lUdrv & Santosll2007h . 

Giant planets, however, seem to be more frequent around 
these stars. The detection limits show that most of the giant plan- 
ets should have been detected in this sample. A frequency of 
2.63% for giant planets around metal-poor stars was calculated 
for stars with at least six measurements, with a sensitivity to pe- 
riods up to 1800 days. If indeed a giant planet was missed in the 
sample, the frequency would be even higher According to sev- 
eral studies (e.g. Santos et al. 2004; Fischer & Valenti 2005), 3% 
is the frequency of giant planets around stars of solar metallic- 
ity. Given the same number, derived in this work, for metal-poor 
stars, this can mean two things: either the planet frequency be- 
comes constant for stars wit h \FelH'\ < 0.0 (for a discuss ion, see 
also e .g. . Santos et al. 2004: lUdrv & Santos 2 007: Johnson et alJ 
l2oToh or previous frequency-models should be higher. 

The metallicities of the stars with discovered giant planets 
all lie above -0.7 dex. Within this metallicity bin, the planet fre- 
quency increases to 4.48;^j 3g%. For a metallicity of -0.55 dex 
and a stellar mass of O.SMq (mean value fo r this sample), previ- 
ous studies report values of fn = 1 .22'^^''^% (iJohnson et alJ2010h 
and fp = 0.14% ( Sousa etal.ll2011b[) .' However, in both cases, a 
powerlaw was fitted over the whole metallicity range (up to 0.6 
and 0.5 dex, respectively). In the low-metallicity end, their fit is 
clearly lower than their observed fraction (~ 4.6% and 3.77%, 
respectively). The value reported here is thus higher than previ- 
ous fits, but consistent within one sigma with the observed frac- 
tion in these previous studies. 

For stars with metallicities lower than -0.7, the frequency 
is lower than 2.36%. In this context, it is worth mention- 
ing that so far, only one (giant) planet has been detected 
around a main-sequence star with metallicity lower than -0.6 
dex dCoc hran et al. 2007). There are some other candidates 
where planets orbit stars with metallicities lower th an -0.6 dex 
dNiedzielski et al.l2009HSetiawan et al.ll2003ll2010l) . These stars 



are all giants, however Furthermore, there is a candidate planet 
detected by imaging that orbits a young main-sequence star with 
a metallicity of -0.64 dex (Chauvin et al. 2005). These planets 
are therefore not relevant for the purpose of this work. 

All the above results are strong evidence that giant planet fre- 
quency is a non-constant function of stellar metallicity. This was 
akeady established for the high-metallicity tail, but this work 
shows that it is also true for lower metallicities. Moreover, the 
frequencies are probably higher than previously thought. A pow- 
erlaw may not be the best function to describe the planet fre- 
quency over the whole metallicity range. As mentioned before, 
this correlation between giant planet frequency and stellar metal- 
licity favors the core-accretion model as the main mechanism for 
giant planet formation. 

With the statistics presented in this work, a metallicity limit 
can be established below which no giant planets can be found 
anymore. According to these statistics, this metallicity limit 
would be about -0.5 — 0.6 dex. This value agrees with the re - 
cent results of the theoretical study of iMordasini et al.l (l2012h . 
These authors looked at correlations between stellar and plane- 
tary properties, based on a synthetic planet population, built by 
the core accretion model. They found that giant planets are not 
formed below -0.5 dex. 

However, more data are still needed to produce better statis- 
tics. Future missions, such as Gaia, will produce better precision 
in the data. For example, in its all-sky global astrometric survey, 
Gaia will probe thousands of nea rby metal-poor stars for gas gi- 
ant planets within 3-4 AU (e.g. ISozzettil|2010l) . thus crucially 
helping to shed light on the planet-metallicity connection. 
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